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Abstract 

The  influence  of  a  set  of  satellite  oscillators  on  the  response  behavior  of  a 
master  oscillator,  to  which  the  set  is  coupled,  is  of  fundamental  significance  to 
structural  acoustics  and  beyond.  The  focus  is  largely  on  the  term  that  modifies 
the  impedance  of  the  master  oscillator  due  to  that  coupling.  The  real  and 
imaginary  parts  in  this  term  are  dubbed  the  induced  reactive  factor  and  the 
induced  loss  factor,  respectively.  For  a  passive  set  of  satellite  oscillators  the 
induced  reactive  factor  may  be  either  positive  (mass-like)  or  negative  (stiffness¬ 
like)  whereas  the  induced  loss  factor  is  invariably  positive  definite.  The  behavior 
of  both  induced  factors  as  functions  of  fi'equency  is  explored,  especially  in  the 
range  of  fi'equency  spanning  the  resonance  frequencies  of  the  satellite  oscillators. 
The  resonance  fi-equency  of  the  master  oscillator  in  isolation  is  placed  centrally 
in  this  frequency  range.  The  phenomenon  of  the  “suppression  of  undulations”,  in 
the  induced  reactive  factor  and  in  the  induced  loss  factor,  as  the  modal  overlap 
parameter  approaches  and  exceeds  unity,  is  demonstrated.  Also  demonstrated 
and  defined  is  the  phenomenon  of  “erosion”  in  these  induced  factors,  again,  as 
the  modal  overlap  parameter  approaches  and  increases  beyond  unity. 
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I.  Introduction 

Harmonic  oscillators  are  simple  structural  elements  that  are  helping  to  model 
structural  complexes.  The  models  are  then  analyzed  in  terms  of  these  coupled 
harmonic  oscillators.  The  result  of  the  analysis  is  interpreted  to  yield  infor¬ 
mation  of  the  sensitivity  of  the  response  behavior  of  the  complex  to  parameters 
that  define  it.  The  sensitivity  may  then  reveal  which  parameters  need  to  be 
modified  in  order  to  achieve,  for  example,  noise  control  goals  for  the  complex. 
The  scheme  of  this  analysis  is  initiated  by  casting  the  structural  complex  in  terms 
of  structural  components.  This  division  of  the  structural  complex  yields  a  matrix 
equation  of  motion  in  the  form 

jf(?  I  jz(x  1  x',t)  =  {jzay  (Xa  I  x'y,t)dx'y),  (I.l) 

and  (0  is  the  temporal  variable.  In  Eqs.  (I.l)  and  (1.2)  the  quantity 
^ay  \^y’0  is  the  element  in  the  impedance  matrix  operator  z{x  \  x',t)  that  is 
due  to  the  (}')th  structural  component  influence  on  the  (a)th  structural  com- 
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ponent.  The  diagonal  elements  are  typically  defined  by  the  self-impedance 
operator  z„{Xf^,t)  and  the  off-diagonal  elements  are  typically  defined  by  the 
coupling  impedance  operator  Zay{Xa\x'y,t)  between  the  {y)th  and  the  {d)th 
structural  components  with  a  there  is  no  coupling  between  these  structural 
components  if  z„y{Xa\Xy,t)  is  identically  equal  to  zero.  Further  in  these 
equations,  v„(x^,0  and  Pecd^a^t)  are,  respectively,  the  response  and  the 
external  drive  for  the  {a)th  structural  component  and  (x^)is  the  vector  that 
spans  the  spatial  dimensionality  of  the  {cc)th  structural  component.  The  rank  of 
the  matrix  equation  of  motion  is  equal  to  {N),  the  number  of  structural 
components  in  the  complex.  Each  structural  component  is  selected  in  a  manner 
that  admits  its  equation  of  motion  to  a  modal  analysis;  i.e.,  the  impedance 
operator  is  selected  to  be  an  eigenoperator  in  the  spatial  domain.  [This 

selection  is  not  necessarily  unique;  there  may  be  a  number  of  choices  for 
selecting  Of  course,  each  choice  may  determine  a  different  set  for 

z^y{x^\x'y,i)  with  a^y  The  modal  analysis  is  characterized  by  eigen¬ 
functions  (mode  shapes),  typically  and  by  eigenvalues,  typically 

Zajit) ,  of  the  eigenoperator  z^{Xc^,t) ;  namely 

(^a » 0  j (^a  ^  ~  ^aj (0  j (^)  •  (f -3) 
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The  mode  shapes  may  be  collated  into  a  mode  shape  vector 

^(^)  “■  {^a(^a)}’  ^a(^a)  ~  (I*^) 

and  the  mode  shapes  are  by  definition,  orthogonal  and  closed 

j.  K 

J  dx^  ^ay(^a)  ^ak^^a)  ~  ^jk  ’  ^^ay(^a)  ~  ~^a)>  (1-5) 

j 

respectively,  where  (^a)  enumerates  the  modes  in  the  (a)th  structural  com¬ 
ponent.  In  Eq.  (I.3)-(I.5)  the  subscripts  (a)  and  (j)  designate  the  (a)th 
structural  component  and  the  (j’)th  mode  in  that  structural  component,  respect¬ 
ively.  The  modes  are  then  employed  to  remove  the  spatial  dependence  from  the 
matrix  equation  of  motion  for  the  complex,  yielding,  in  its  place,  a  modal  matrix 
equation  of  motion  in  the  form 

f(Ov(0  =  Peit)  ,  (1.6) 

where  z(t)  is  the  modal  impedance  matrix  operator,  y(r)  is  the  modal  response 

vector  and  Pg(t)  is  the  modal  external  drive  vector.  In  some  details 

N  Ny 

|(0  =  (|aa  fa/  ~  fa/  (i  ~  fa/))»  |aa  “  {.^aj  ^jk)’  ^aj  ~  ^^^yqaj » 

r  q 

^ajai  =  0  for  i;  Z^a  =  i^aJrq)’ 

^ajyq  ~  J^a  ^aj^af  ■^a/('^a  ^^y  ^/^(^/)> 

V  =  {fa } »  fa  =  {Vay  } ;  (x^ ,  t)  <l>^j  ) ,  (I.7b) 

Pe  ~  ^Pecc }  >  P^(^  ~  ^Pecc  j^’  Pea  j  ~  J^^a  Pea  (■^a  ’  ^a  j  (^a  )  > 


(I.7c) 


6 


where  occasionally,  the  dependence  of  quantities  on  the  temporal  variable  is 
obvious  and,  as  such,  is  suppressed.  The  particular  division  of  the  structure  into 
components,  each  of  which  admits  to  a  modal  decomposition,  defines  the  nature 
of  the  couplings  among  the  components.  Clearly,  the  orthogonality  of  the  modes 
in  a  structural  component  assures  that  modes  that  belong  to  that  component  are 
not  coupled  to  each  other;  i.e.,  z^a » in  Eq.  (I.7a),  is  a  diagonal  matrix  operator. 
Moreover,  the  modal  impedance  operator  is  equivalently  that  of  a  harmonic 
oscillator.  The  coupling  between  two  structural  components  is  then  defined  in 
terms  of  coupling  between  a  mode  of  one  component  to  that  in  another; 
Zayd-^a/)  accounts  for  off-diagonal  elements  only.  For  example,  the  off- 
diagonal  element  Zaj/g  pertains  to  the  coupling  element  between  the  (q)th 
mode  in  the  structural  component  and  the  (J)th  mode  in  the  (a)th 

structural  component  with  .  The  rank  of  the  modal  equation  of  motion  is 

N 

equal  to  where  (iV^)  is  the  number  of  viable  modes  in  the  (a)th 

a 

structural  component.  Naturally,  the  rank  of  Eq.  (1.6)  usually  far  exceeds  that  of 
Eq.  (I.l);  that  increase  in  rank  is  the  penalty  that  is  paid  for  the  removal  of  the 


spatial  dependence  from  the  latter  equation.  The  penalty  is  partially  relieved  by 
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the  fact  that  the  modal  matrix  impedance  operator  in  Eq.  (1.6)  harbors  elements 
that  are  much  simpler  than  those  usually  characterizing  the  matrix  impedance 
operator  in  Eq.  (I.l),  However,  the  increase  in  rank  is  not  the  only  penalty.  If 
spatial  information  in  the  response  of  the  structure  is  significant,  to  reconstitute 
this  information  one  requires,  in  addition  to  the  derivation  of  y(0,  the  construct¬ 
ion  of  the  scalar  product 

’ 0  ~  liz  (0  *  (•^a  ) ’ 

where  y(x,/),  y(0  and  y^if)  are  defined  Eqs.  (1.2),  (1.4),  (1.6)  and 

(I.7b),  respectively,  and  the  superscript  (T)  indicates  the  transpose  of  the 
quantity.  The  availability  and  storage  of  the  mode  shape  vector  ^(x),  as  stated  in 
Eq.  (1.4),  is  a  tall  order  that  can  seldom  be  satisfactorily  met.  Whether  the 
penalty  for  circumventing  Eq.  (I.l)  via  Eqs.  (1.6)  and  (1.8)  is  advantageous,  is 
difficult  to  assess  without  clear  statements  of  goals  and  requirements.  Often  such 
statements  can  be  made  only  when  the  analysis  is  completed;  a  catch  (22),  to  be 
sure.  Excluded  from  these  arguments  are  situations  in  which  the  analyzed 
complex  ideally  lacks  a  priori  spatial  dependence  and  yet  the  constituent 
components  can  be  represented  by  sets  of  harmonic  oscillators.  In  this  case  Eq. 


(1.6)  is  the  whole  shebang. 
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In  this  thesis  a  modest  part,  in  this  grand  scheme  of  providing  descriptions  to 
the  behavior  of  structural  complexes,  is  pursued.  The  modal  matrix  equation  of 
motion  is  assumed  given.  The  aim  is  to  seek  answers  to  questions  such  as:  How 
is  one  to  draw  relevant  information  that  lies  within  the  modal  equation  of 
motion?  And  how  is  one  to  interpret  this  information  in  seeking  to  beneficially 
modify  the  parameters  that  specify  the  complex  that  is  at  the  core  of  this  equation 
of  motion?  Moreover,  the  thesis  deals  with  a  highly  simplified  complex 
composed  of  a  single  master  oscillator  coupled  to  a  set  of  satellite  oscillators. 
The  satellite  oscillators  in  the  set  are  uncoupled  to  each  other.  Several  papers 
and  theses  are  found  in  the  open  literature  with  a  largely  similar  aim  in  sight, 
employing  basically  a  similar  complex  [1-11].  What  is  then  new  in  the  present 
effort?  For  one,  most  of  these  cited  papers  deal  with  satellite  oscillators  that  are 
essentially  sprung-masses  as  depicted  in  Fig.  la.  The  coupling  in  this  case  are 
stiffiiess  control  and  they  are  strong.  In  the  present  thesis  other  forms  of 
couplings  are  also  introduced,  in  particular,  mass  and  gyroscopically  controlled 
couplings  [12,  13].  In  addition,  the  coupling  strengths  are  moderated  so  that  the 
influences  of  moderate  and  weak  couplings  may  be  defined  and  analyzed.  The 
modal  equation  of  motion  is  stated  for  the  complex  as  depicted  in  Fig.  lb.  Mass, 
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stiffiiess  and  gyroscopic  coupling  parameters  are  defined  and  employed  to 
establish  the  equation  of  motion  for  a  complex  that  harbors  such  forms  of 
couplings. 

Assuming  that  the  complex  is  stationary  in  the  temporal  domain,  it  is 
convenient  to  express  the  modal  matrix  equation  of  motion,  hereto  defined  in  the 
time  domain,  also  in  the  frequency  domain.  In  the  frequency  domain  Eqs.  (1.6) 
and  (1.7)  take  the  forms 


(1.9) 

^(fi>)  =  (_^aa(^)^ay  ~  ?aa  ~  (.^aj  ^jk)’ 

N 

^ajai  =0  for  7^1 5 

r  q 

/(() = /(/), 

(1. 10a) 

v  =  {Z.};  e,  = 

(1. 10b) 

(1. 10c) 

where  the  Fourier  function  / (t)  is  defined 

/(O  = 

(1. 11) 

and,  again,  occasionally  the  dependence  of  quantities  on  the  frequency  variable  is 
obvious  and,  as  such,  is  suppressed.  The  choice  of  the  frequency  domain  versus 
the  temporal  domain  is,  with  caution,  largely  a  matter  of  convenience  and 
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preference;  again  with  caution,  the  two  domains  are  merely  related  by  a  Fourier 
transformation  [14],  One  recalls  that  the  modal  matrix  equation  of  motion  is 
defined  in  terms  of  a  matrix  impedance  operator  that  describes  the  behavior  of  a 
complex  composed  of  coupled  harmonic  oscillators.  In  the  model  of  the 
complex  here  considered,  a  master  oscillator  is  coupled  to  a  set  of  satellite 
oscillators;  the  satellite  oscillators  are  not  coupled  to  one  another,  [cf  Fig  lb.] 
The  modification  in  the  matrix  element  that  describes  the  impedance  operator  of 
the  master  oscillator,  induced  by  the  coupling  to  the  set  of  satellite  oscillators,  is 
of  particular  interest.  This  modification  is  characterized  by  a  reactive  part  S(y) 
and  by  a  real  part  where  (y)  is  the  normalized  frequency  variable; 

y  =  {o)/o)o),  with  {cOg)  the  resonance  frequency  of  the  master  oscillator  in  iso¬ 
lation.  A  positive  value  for  the  induced  reactive  factor  S{y)  contributes  a  mass 
control  term,  a  negative  value  contributes  a  stiffness  control  term.  For  the 
passive  complex  here  analyzed,  the  induced  loss  factor  ri^{y)  is  positive  definite. 
The  natures  of  S{y)  and  of  7^(y)and  their  modifications  to  variations  in  the 
parameters  that  specify  the  complex  are  recorded  and  interpreted.  Some  of  these 


variations  are  locally  imposed  by  pseudo-statistical  means  e.g.,  the  resonance 
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frequency  and  the  modal  overlap  parameters  that  are  assigned  to  each  individual 
satellite  oscillator  are  picked  pseudo-statistically,  either  separately  or  in  unison. 

A  companion  report  that  deals  only  with  the  induced  loss  factor  ;75(>^)has 
been  recently  issued;  in  that  report  the  induced  reactive  factor  S(y)  is  not  invest¬ 
igated  [15].  Indeed,  the  present  report  is  largely  issued  in  order  to  rectify  this 
omission.  In  the  present  report  it  is  noted  that  the  induced  reactive  factor  S(y) 
may  be  of  significance  only  when  the  coupling  strengths  are  strong.  In  these 
cases  at  the  edges  of  the  frequency  range  spanned  by  the  resonance  frequency 
distribution  of  the  satellite  oscillators,  the  absolute  value  of  the  induced  reactive 
factor  S(y)  may  approach  and  exceed  unity.  On  the  other  hand,  even  at  the 
edges  of  this  frequency  range,  for  moderate  coupling  strengths  the  values  of 
1 5(3;)  I  are  insignificant  and  for  weak  coupling  strengths  the  values  of  |  S(y)  |  are 
negligible  compared  with  unity. 

A  loss  factor  is  basically  an  energetic  quantity,  and,  therefore,  one  may 
inquire  whether  the  induced  loss  factor  7^(7),  as  just  determined  via  the  linear 
impedance  analysis  (LIA)  may  be  estimated  via  an  energy  analysis  (EA).  In  this 
vein  in  another  companion  report,  the  analysis  of  the  complex  performed  herein 
and  in  Reference  15  is  repeated  [16].  In  Reference  16,  however,  instead  of  using 
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a  linear  impedance  analysis:  e.g.,  as  stated  in  Eq.  (1.9),  an  alternative  approach  is 
used.  In  this  approach  the  stored  energy  E^iy)  in  the  master  oscillator,  the 
stored  energy  Ej.{y)  in  the  {r)th  satellite  oscillator  and  the  stored  energy  E^^^y) 
in  the  coupling  of  this  satellite  oscillator  to  the  master  oscillator  are  estimated  in 
terms  of  the  sum  of  the  kinetic  and  potential  stored  energies  in  each  category. 
Distinguishing  between  the  kinetic  and  potential  stored  energies  renders  it 
possible  to  determine  the  Lagrangians  for  the  complex.  Applying  the  Lagrangian 
equations,  the  linear  equation  of  motion  are  yielded  [12,  16].  This  procedure 
provides  a  check  that  the  equation  of  motion  in  the  linear  impedance  analysis 
(LIA)  are  compatible  with  the  equation  of  motion  in  the  energy  analysis  (EA). 
This  check  is  successfully  performed  [12,  16].  In  terms  of  the  stored  energies, 
^r(T)  and  E^riy),  the  dissipated  power  no(>;)  in  the  master  oscillator, 
the  dissipated  power  n^(;/)in  the  {r)th  satellite  oscillator  and  the  dissipated 
power  Tl^j.(y)  in  the  coupling  between  these  two  oscillators  are  estimated, 
respectively,  using  the  loss  factors  ,  rj^  and  rj^r  that  are  assigned  to  the  master 
oscillator,  to  the  {r)th  satellite  oscillator  and  to  the  coupling  between  them, 
respectively.  Also  estimated  is  the  external  input  power  Hgiy)  into  the  master 
oscillator  by  an  external  drive.  Under  the  prevailing  assumptions  that  only  the 
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master  oscillator  is  driven,  Tlg(y)  is  the  total  power  input  into  the  complex  from 
an  external  drive.  Invoking  the  conservation  of  energy  (or  the  balance  of  power) 
one  may  equate 

=  n„(>;)+n^,(3;);  n^,(y)  =  n,(;;)  +  n^,(y); 

N 

no5(>^)=Enor(3^).  (I.12a) 

r=l 

where  (N)  is  the  number  of  satellite  oscillators  in  the  complex.  From  Eq.  (1. 12a) 
the  loss  factor  77^(3^)  of  the  master  oscillator  and  the  induced  loss  factor  Tjg{y) 
may  be  estimated;  namely 

(y^/2)[n^(y)/coE^K(y)]  =  ^o(y)  ^  (I.12b) 

{y^ l2)[Uo,iy)la)EgKiy)}  =  vtiy)  ,  (I.12c) 

where  E^j^  (3/)  is  the  kinetic  energy  stored  in  the  master  oscillator,  as  a  function 
of  (y ) .  The  superscript  (e)  is  introduced  to  distinguish  between  the  loss  factor 
77o  of  the  master  oscillator  and  the  induced  loss  factor  (y) ,  as  assigned  and 
determined  for  the  master  oscillator,  respectively,  via  the  linear  impedance 
analysis  (LI A),  and  the  loss  factor  70(37)  and  the  induced  loss  factor  rj^iy),  as 
determined  via  the  energy  analysis  (EA),  respectively.  The  equality  of  the  loss 
factors  tj^  and  7® (7)  is  obvious  and  is  stated  as  such  in  Eq.  (1. 12b).  On  the 
other  hand,  the  comparison  between  the  two  induced  loss  factors;  7^  (37)  and 
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rj^iy) ,  is  computationally  performed  and  it  is  found  that,  by  and  large,  they  are 
identical  [16]. 

Finally,  and  in  addition,  in  Reference  16  the  relationship  of  (EA)  to  the 
statistical  energy  analysis  (SEA)  is  established.  In  particular,  it  is  found  that 
(SEA)  may  underestimate  the  ratio  H;^(y)  of  the  energy  EQ^{y)  stored  in  the 
satellite  oscillators  and  in  the  couplings  to  the  energy  E^^y)  stored  in  the  master 
oscillator;  namely 

N 

HoW  =  [£„W/4W];  EM  =  Hem-^ 

r=\ 

Eoriy)  - Eriy)  +  E^r{y)\  E^{y)  =  E„i^(,y)  +  Eop{y),  (I.13a) 

where  {N),  again,  is  the  number  of  satellite  oscillators  and  E^piy)  is  the  potent¬ 
ial  energy  stored  in  the  master  oscillator,  as  functions  of  (y) .  [cf  Eq.  (I.  12).] 
The  ratio  H^(;j;)  is  a  measure  of  the  global  coupling  strength  in  the  milieu  of 
(SEA);  it  is  related  to  the  modal  coupling  strength  (y)  in  the  form 

K(y)  =  K(t)/«.(t)]  Co(y),  (i.i3b) 

where  %(y)  and  Kj(3^)  are  the  modal  densities  in  the  master  oscillator  and  in  the 
set  of  satellite  oscillators,  respectively.  (A  modal  density  is  the  average  number 
of  modes  per  unit  frequency  [12].)  A  tenet  in  (SEA)  is  that  the  modal  coupling 
strength  (t)  can  never  exceed  unity.  It  transpires,  nonetheless,  that  for  the 
stronger  coupling  strengths,  the  modal  coupling  strength  estimated  via 

(EA)  may  exceed  unity  when  the  typical  modal  overlap  parameter  (b),  assigned 
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to  the  satellite  oscillators,  falls  below  a  threshold  value.  Thus  a  discrepancy 
between  the  modal  coupling  strengths,  as  estimated  via  (EA)  and  via  (SEA),  may 
ominously  emerge.  However,  for  modal  overlap  parameters  that  exceed  the 
threshold  value,  the  validation  of  (SEA)  is  not  challenged  by  (EA)  [16]. 
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11.  Elaborated  Coupling  Forms 

In  a  companion  paper,  designated  as  Reference  1,  the  complex  is  specified  in 
terms  of  the  impedance  Zq(<»)  of  the  isolated  master  oscillator,  the  impedance 
(/ft)w^)of  the  (r)th  isolated  satellite  “oscillator”  and  the  coupling  impedance 
{kj.  I  ico);  the  coupling  is  between  the  master  oscillator  and  the  satellite  oscillator. 
This  complex  is  sketched  in  Fig. la  [1-11].  In  this  complex,  the  satellite 
oscillators  remain  uncoupled  to  each  other.  The  purpose  in  the  present  paper  is 
to  introduce  a  significant  extension  in  scope.  In  this  extension  the  isolated 
master  oscillator  remains  the  same.  An  isolated  satellite  oscillator,  the  (r)  th,  is 
specified  by  the  impedance  Z^(iy)  that  may  consist  of  both  a  mass  and  a 
stiffiiess  control  term.  The  coupling  of  this  satellite  oscillator  to  the  master 
oscillator  is  specified  by  the  coupling  impedance  and  a  gyroscopic 

coupling  coefficient  (G^.)  [12,13].  The  coupling  impedance  Z^f.{(a)  consists  of  a 
mass  and  a  stiffiiess  control  term.  A  complex  of  this  type  is  sketched  in  Fig.  lb. 


The  mass  and  the  stiffness  control  terms  in  this  complex  are  related  in  the  forms 
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K  =K,(l+i?j,); 

{KJM)  =(ol. 

(la) 

K  =^or(l  +  ^7r); 

iKrl^r) 

(lb) 

^cr  ~  kocr^'^^Vcr^ 

’  ^^ocr  ^  ’ 

(Ic) 

where  the  pairs  {M,K),  and  {mj.,kcr}  are,  respectively,  in  reference 

to  the  master  oscillator,  to  the  ir)th  satellite  oscillator  and  to  the  coupling 
between  them.  The  parameters  (tJo),  (77^),  and  {rjcr)  are  the  corresponding 
stiffiiess  control  loss  factors,  respectively.  To  complete  the  definition  of  the 
coupling,  the  mass  (m£.^)and  the  gyroscopic  coefficient  (G^.)  need  to  be 
specified  [12,  13].  These  are  specified  through  the  coupling  parameters  which 
are  defined  in  the  normalized  forms 

^cr  =  Kr  l^r)  ^  fr  =  [<^r  > ^r)]  >  (Id) 

respectively.  With  the  assistance  of  Fig.  lb,  the  linear  equations  of  motion  of  the 

master  oscillator  in  situ  and  of  a  typical  satellite  oscillator  in  situ  are  derived 
R 

Z°(<»)F<,(©)+y;z„(®)  +  [Z- (ffl)  -  (2a) 

1 

[Z,  (co)+Z,,(a))]V,(co)  +  [Z- (fi>)+G,  ]F^(o)  =  0,  (2b) 

respectively,  where  Vg(a))  and  (to)  are  the  responses  of  the  mass  (M)  of  the 
master  oscillator  and  the  mass  (m^)  of  the  (r)th  satellite  oscillator,  respectively. 
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(R)  is  the  number  of  satellite  oscillators  that  are  coupled  to  the  master  oscillator, 
is  the  drive  that  is  applied  externally  to  the  master  oscillator;  the  satellite 
oscillators  are  not  driven  externally,  Z^(o))  is  the  impedance  of  the  master 
oscillator  in  isolation 

Zo(o>)  =  (i(oM)[l-{y)~^(\+iTj^)];  y  =  (a)to)„),  (3a) 

is  the  impedance  of  the  (r)th  satellite  oscillator  in  isolation 
=  +  z^=(x^/y),  (3b) 

Z^^{a>)  is  the  impedance  of  the  coupling  between  the  master  oscillator  and  the 
(r)th  satellite  oscillator 
Zcr  ((o)  =  (i(om,)  -  {z^,  f  (1 + /  )] ; 

^cr=Kr/w^);  Xcr={0)crlo>o)\  Z^,={Xcrly),  (3c) 

Z^^(<u)  is  an  impedance  that  is  related  to  the  coupling  impedance  Z^^(fi)); 
namely 

^c"r(<w)  =  )[^cr  +(^cr)^  (1  +  ^  )] »  (3d) 

and  {fhcr),{Xcr\{ricj.\  and  (G^)are  defined  in  Eqs.  (1)  and  (2)  [12].  The 
superscript  (o)  is  reserved  to  designate  quantities  that  pertain  either  to  the  master 
oscillator  in  isolation  or  to  satellite  oscillators  under  certain  and  definitive 
impositions.  Thus,  for  example,  from  Eqs.  (2a)  and  (3a),  one  may  state 


19 


Z^(a?)Ff(a?)  =  P,(a?),  (4) 

where  V^(a))is  the  response  of  the  master  oscillator  in  isolation.  By  a 

straightforward  algebraic  manipulation  of  Eq.  (2)  one  derives 

Zo(£0)  V^(ci>)  =  ;  F^Ccd)  =  B^(ci?)F^(a?) ,  (5) 

where 

R 

Z„(®)  =  z“(<i))  +  '£l{Z,(<»)Z„Ca>)}+(Q^f][Zr(<»)+Zcr(‘»)r',  (6a) 

1 

^r(^)  =  [Zcr(^)  +  (6b) 

(Qcrf=4m^,k^+(G,f,  (6c) 

and  the  quantities  Zq^co),  Z^(a)),  Z^^(6))and  Z^(<»)are  stated  explicitly  in  Eq. 
(3).  Indeed,  from  Eqs.  (1),  (3)  and  (6)  one  obtains 

Z^((o)  =  Z^(y)  =  (i(oM)[l-(y)-^{[l-S(y)]+i[;/^  +  7/^(y)]}],  (7a) 

5^(<u)  =  -[m^^  +  (z^rf  (1  -airier)  -  iigr  ly)]  • 

[(1  +  frier)- i^rr  f  (1  +  i Irr )]“* .  (7b) 

where 

(Xrr)  (1  +  ir/rr)  =  (Xr)  (l  +  ZT/cr)  (^cr) 

R 

[5(3;)  -i?J, (7)]  =  (y)^  {iHr  { [1  -  {Zr  f  (1  yirir)]* 

1 

[nier  -  {.Zer  f  (1  +  i  TJer  )]  “  (^cr  !  yf)  ^ 

[(1  ^fher)-  {Zrr  f  (1  +  i Tlrr )r^ } , 


(8b) 
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and,  again 

^rr  ~  (^rr  ^3^)  s  ~  ^3^)  j  ^cr  ~  i^cr  ^3^)  > 

(.^cr^yf  =  (^crf  (1  +  i  Vcr)  +  (^r  /  ;  ^cr  =  iQcr  /K  '"r  )]•  (8c) 

One  notes  that  the  compound  coupling  parameter  is  a  functional  of  the 
mass  and  gyroscopic  coupling  parameters  (w^^)and  (g^.),  respectively.  These 
coupling  parameters  are  defined  in  Eq.  (Id).  One  also  notes,  with  satisfaction, 
that  the  dependence  of  the  terms  in  the  sum  on  the  gyroscopic  coupling 
parameter  (g>)  is  quadratic  so  that  the  sign  assigned  to  the  gyroscopic  coefficient 
{Gf.)  plays  no  role  in  the  influence  of  the  individual  satellite  oscillators  on  the 
impedance  of  the  master  oscillator.  The  gyroscopic  coupling  is  in  quadrature  to 
both,  the  mass  and  the  stiffness  control  couplings. 

Examination  of  Eqs.  (7)  and  (8)  shows  that  the  normalized  impedance  that  the 
satellite  oscillators  collectively  induce  on  the  master  oscillator  may  be  cast  in 
term  of  the  two-vector  {^(y),  75(7)}  (/?),  which  is  a  function  of  (R),  as 
indicated.  The  evaluation  of  this  two-vector,  however,  is  predicated  on  explicitly 
specifying  the  two-vector  7^^ }(/?),  its  two  supplemental  components 
{x^,  7^}(i2)  and  {x^.;.,  7c;.}(i?)and  finally  assigning  the  compound  coupling 
(jqcr)-  The  two-vector  {x^^,  7^^}  (i?)is  designed,  for  the  sake  of  convenience,  to 
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stay  fixed  with  respect  to  variations  in  the  index  (r) .  In  this  design,  the  springs, 
that  are  placed  on  either  side  of  the  mass  of  a  satellite  oscillator,  are  set  to  be 
similar.  With  this  design,  Eq.  (8b)  is  evaluated.  In  this  paper  the  induced 
reactive  factor  iS'(y)  and  the  induced  loss  factor  rj^iy)  are  considered.  These  two 
quantities  are  examined  for  a  variety  of  coupling  forms  and  coupling  strengths  as 
well  as  for  a  number  of  values  of  the  modal  overlap  parameters  (3^)  and 
associated  with  the  loss  factors  7^  and  7^^ ,  respectively  [1  ].  By  and  large,  (Z>^)  is 
set  equal  to  (Z>c^)and  they  are  equated  to  {h)-,  b  =  br  =  b^r.  The  simplifying 
equalities  are  imposed  without  considerable  loss  in  generality.  Using  Eq.  (8b), 
the  exact  evaluations  of  S(y)  and  of  ?Js(y)  are  executed  for  three  values  of  (Z>); 
b  =  (0.1),  (2.0)  and  (10) .  In  these  evaluations  one  finds  that  S{y)  and  rig{y)  are 
functions  of  (b).  On  the  other  hand,  the  first  order  approximation  of  iS'(7)  and  of 
Vsiy)  5  designated  hy Sj(y)  and  by  7/ (y),  respectively,  that  are  derived  from  the 
replacement  of  the  summation  in  Eq.  (8b)  by  an  integration,  are  found  to  be 
independent  of  (b).  The  evaluations  are  graphically  displayed;  three  exact  eval¬ 
uations  for  the  three  values  of  (b)  and  the  corresponding  first  order 
approximation  of  5(3;)  and  7,(7)  are  superimposed  in  each  of  the  displays.  The 
curves  representing  the  first  order  approximations  of  S(y)  and  of  tj^iy) ,  in  these 
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displays,  are  dubbed  (FOA).  The  comparison  between  these  four  evaluations  in 
each  display  can  thus  be  made  at  a  glance,  assisting  greatly  in  the  interpretation 
of  the  data  that  is  computed  and  displayed.  In  the  frequency  range  of  concern, 
when  (b)  is  small  compared  with  unity  the  values  of  S(y)  and  the  levels  of 
7]s(y)  in  the  exact  computations,  undulate,  as  functions  of  (y).  (The  frequency 
range  of  concern  spans  the  resonance  frequencies  of  the  satellite  oscillators.) 
The  excursions  in  the  undulations  increase  with  decrease  in  (b).  It  transpires  that 
the  mean-value  averaging  of  the  undulated  values  of  S(y)  and  levels  of  tj^iy) 
when  (b)  is  small  compared  with  unity;  b«l,  are,  respectively,  coincident  with 
curves  that  are  independent  of  (b)  [15-17].  (Values  are  arithmatically  averaged, 
levels  are  geometrically  averaged.)  These  curves  for  S(y)  are  designated  iS'^(>’) 
and  for  77^(7)  are  designated  Jj^(y).  In  the  frequency  range  spanned  by  the 
resonance  frequency  distribution  of  the  satellite  oscillators  -  -  the  frequency 
range  of  concern  -  -  S^(y)  is  positive;  i.e.,  mass  control,  for  y<^  and  is 
negative;  i.e.,  stiffiiess  control,  for  y>  1 .  Moreover,  in  that  frequency  range  of 
concern,  S^iy)  is  largely  a  “straight-line”  intersecting  the  neutral  line  either  at  or 
in  the  vicinity  of  >>=1 .  On  the  other  hand,  in  that  same  frequency  range  rj^iy)  is 


also  a  “straight-line”;  that  line,  except  when  a  gyroscopic  coupling  is  present,  lies 
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largely  horizontal  at  a  given  level.  As  (b)  approaches  and  increases  above  unity 
the  undulations  in  the  exact  values  of  S{y)  and  in  the  exact  levels  of  are 
suppressed  and,  in  the  frequency  range  just  defined,  the  phenomenon  of  erosion 
commences  and  increases  with  increasing  {b) .  In  the  case  of  the  induced 
reactive  factor,  erosion  is  manifested  by  values  that  progressively  converge  from 
.5^(3;)  onto  the  first  order  approximation  to  S(y),  as  (b)  increases.  The  first 
order  approximation  to  S(y);  namely  Sj(y),  is  equal  to  zero;  Sj(y)=0  and  in 
general,  S^iy)  ^  Sj(y).  On  the  other  hand,  in  the  case  of  the  induced  loss 
factor,  erosion  is  manifested  by  levels  that  progressively  decrease,  from  the  first 
order  approximation  of  7^(y),  as  (6)  increases  [6,  7].  The  first  order  approx¬ 
imation  of  TJsiy) ;  namely  7jj{y) ,  is  coincident  with  Tj^iy);  rj^iy)  =  7/(y) . 
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III.  Resonance  Frequency  Distribution  of  and  Assignment  of  Individual 
Loss  Factors  to  the  Satellite  Oscillators 

It  has  been  argued  that  the  coupling  between  two  oscillators  may  cause  a  shift 
in  the  resonance  fi'equency  that  each  has  in  isolation.  Indeed,  it  was  suggested 
that  this  shift  may  be  used  to  determine  the  “coupling  strength”  [18,  19].  One 
may  then  question  as  to  what  exactly  are  these  shifts  and  are  they  significant. 
The  shifts,  as  such,  and  the  implications  that  they  may  harbor  are  not  addressed 
in  the  present  paper;  here  these  shifts  are  overridden  by  design.  The  design 
intends  to  derive  a  suitable  resonance  frequency  distribution  and  a  proper 
assignment  of  individual  loss  factors  to  the  satellite  oscillators.  The  design  is 
expressed,  then,  in  terms  of  the  two-vector  where  (x^;.)is  the 

normalized  resonance  frequency  and  {rij.^)\s  the  loss  factor  associated  with  the 
ir)th  satellite  oscillator.  Examination  of  Eqs.  (7)  and  (8)  shows  that  the  nor¬ 
malized  resonance  frequencies  of  the  satellite  oscillators,  in  situ,  are  ascertained 
by  satisfying  the  equality 

{^  +  Kr)iyf  =i^rrf‘^  (1  +  +  {x^rf .  (9a) 


As  in  Reference  1,  here  too  (x^^)is  assigned  a  priori  with  equal  numbers  of 
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resonance  jfrequencies  on  either  side  of  the  resonance  frequency  of  the 
master  oscillator  and  the  distribution  is  aligned  in  ascending  order;  namely 
x,,<Xqq-,  q  =  {r  +  \)\  l<r<(i?-l).  (9b) 

A  scheme  that  simultaneously  satisfies  Eqs.  (9a)  and  (9b)  and  the  mid-point  just 
imposed,  demands  a  supplemental  condition  of  design  between  the  two  spring 
stiffiiesses  that  support  the  mass  of  a  satellite  oscillator,  [cf.  Fig  lb.]  This 
condition  of  design  requires  that 

{^rrf  =  (1  +  ^cr)  ^1“  F  =  (1/2)  ,  (9c) 

where,  in  Eq.  (9c),  (F)  may  be  allowed  a  continuous  connotation  [15,  16].  Con¬ 
sulting  Reference  1  and  using  Eq.  (9a),  one  may  express  the  loss  factors  that  are 
associated  with  the  {r)th  satellite  oscillator  in  the  forms 
nr  =  \.M2(x,f\d(x,fidrv,  =[b„n(x„f}{d(.x„f  /dry, 

^rr  ~  [(■^rr)  1  {(^r)  (^cr)  Vcr'^’  (1®) 

where  (/>^)  and  (/>^^)  are  the  modal  overlap  parameters  assigned  to  the  back 
spring  and  the  fore  spring.  These  two  springs  support  the  mass  (m^). 
[cf.  Fig.  lb.]  For  the  sake  of  simplicity  and  computational  advantage,  the 
stiffiiess  on  either  side  of  the  mass  of  a  satellite  oscillator  are  assumed  to  be 


“similar”  in  the  sense 
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(x,)  =  (x,,)  =  (x^);  (x,,)  =  {a,  +  a,,f^ix^),  (11) 

where  (a^)  and  (a^^)are  dubbed  the  spring  factors.  In  keeping  with  the 
definitions  of  (m^r)  (.Sr)’  the  mass  coupling  parameter  and  the  gyroscopic 

coupling  parameter,  the  spring  factor  (a^^)  may  be  designated  the  stiffness 
coupling  parameter.  It  is  conducive  to  specify  (x®)  in  a  form  that  is  compatible 
with  one  of  the  forms  introduced  in  Reference  1;  namely 
x;=[l+{l-2F}K^)^‘'^  (12) 

where 

F  =  riR  +  ir^;  R  =  R(R  +  ir^;  r(R)  =  [r/(2R)];  r<h  (13) 

and  (r)  may  be  discrete;  1  <  r  <  ,  or  continuous  (f)  <  r  <  (i?  +  s)  with  s  <1. 

It  should  be  appreciated,  however,  that  although  in  this  report  Eq.  (12)  is  cast  in 
stone,  other  forms  for  (x°)may  be  readily  introduced  and  similarly  manipulated. 
Also  and  similarly,  as  introduced  in  Reference  1,  the  normalized  mass  (m^)of 
the  (r)th  satellite  oscillator  is  assumed  to  be  independent  of  (r)  and  to  be  of  the 
form 

R 

m,=im,/M)  =  (M,/MXR)-';  M,  (14) 

1 

With  the  intended  exceptions  of  the  last  section  in  this  report  it  is  convenient, 
without  a  great  loss  in  generality,  to  assume  that  the  spring  factors  (a^)  and 
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(a^r)’  the  coupling  parameters  (g^)  and  and  the  modal  overlap  parameters 
(br)  and  (b^^)  are  to  be  independent  of  (r) ;  namely 

ar=a;  a^r^a^;  gr=g;  mcr=^c>  b,=bcr=b,  (15) 

and  it  is  observed,  as  already  intimated,  that  the  modal  overlap  parameters  (b^) 
and  are  set  equal  to  (b).  [cf.  Appendix  A.]  From  Eqs.  (9c),  (10),  (11)  and 


(12),  one  then  finds 

(« +  «c)  =  (l  +  ^c)>  (9d) 

=  Icr  =  Vrr  =  WY  =  /^)  {r{R)ix°f  ],  (16) 

where 

6  =[(;rZ>)  (/?  +  !)“*].  (17) 

Under  this  imposition,  the  two-vector  {(x^^),  {iJrr)}  {R)  assumes  the  simple  form 
inrr)}  (R)  =  {(1  + (x®),  7(F)},  (18a) 

and,  if  further,  the  mass  coupling  parameter  (m^)  is  negligible;  (m^)  «  1 ,  then 

{(4).  (7rr)}  (R)  =  mh  (18b) 


where  the  superscript  (o)  in  (x^^)  and  (7^^)  recognizes  that  (x°^)  and  (7®^)  are 
restricted  to  specific  impositions,  [cf  Eq.  (4).]  With  R=  27  and  =  0,  Eq. 
(18)  is  evaluated  and  depicted  in  Fig.  2a,  as  a  function  of  (F) ;  Fig.  2a.  1  depicts 
(x^^)  and  Fig.  2a.2  depicts  (7^^.) .  In  Fig.  2a  the  modal  overlap  parameter  (b)  is 
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increased  from  (0.1)  to  (2.0)  and  then  onto  (10).  To  these  changes  (77^^)  increases 
by  a  factor  of  (20)  and  then  by  a  factor  of  (10^),  whereas,  to  these  changes  in 
(6),  (x^^)  remains  intact.  With  R  =  7andmc=0,  Eq.  (18)  is  evaluated  and 
depicted,  in  the  format  of  Fig.  2a,  in  Fig.  2b.  On  the  other  hand,  with  R  =  27, 
but  with  Wc=0.75,  where,  for  example,  in  addition  a  =  1.75,  =  0,  Eq. 

(18)  is  evaluated  and  depicted,  in  the  format  of  Fig.  2a,  in  Fig.  2c.  It  is  noted  that 
for  a  continuous  r ,  except  for  obvious  end  conditions.  Figs.  2a.  1  and  2b.  1  are 
identical.  However,  Figs.  2a.2  and  2b.2  are  not  identical.  On  the  other  hand,  it  is 
noted  that  Figs.  2a.2  and  2c.2  are  identical,  however,  Figs.  2a.  1  and  2c.  1  are  not 
identical.  The  commonalties  and  disparities  among  Figs.  2a,  2b  and  2c,  are  as 
expected.  [It  is  also  noted,  in  passing,  that  {x^r)  and  (t/^.;.),  stated  in  Eq.  (18)  are 
independent  of  the  gyroscopic  coupling  parameter  (g;.).] 

It  remains  then  to  use  the  two- vector  stated  in  Eq.  (18)  to  evaluate  the  induced 
reactive  factor  S{y)  and  the  induced  loss  factor  T]g{y) .  Indeed,  using  Eq.  (18)  in 
Eq.  (8b),  one  derives  the  more  explicit  expression  for  S{y)  and  for  7^  (7)  in  the 


form 
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R 

[S(y)-iVs(y)]  =  (yf  W,-a,{z°f{Uirj{r)}] 

1 

-Afh^  {z°f  {\^iW)}-{glyf  }  •  [(l+m^)-(a+ac)(^r  ; 

z^={x°ly),  (8d) 

where  (w^),  (g),  (;;),  {a),  (a^),  {x°),  (w^)  and  7(r )  are  stated  in  Eqs.  (Id),  (3a), 
(11),  (12),  (14)  and  (16),  respectively.  The  computations  of  S(y)  and  of  rj^iy) 
are  largely  carried  out  assigning  the  standard  values 

(M,  I M)  =  10"^ ,  6  =  (0.1)  and  R  =  27  ,  (19) 

where  {M^)  is  stated  in  Eq.  (14),  (b)  is  the  modal  overlap  parameter  and  (R)  is 
the  number  of  satellite  oscillators  in  the  set.  When  these  standard  assignments 
are  deviated  from,  specific  mentions  are  to  be  rendered,  notwithstanding  that,  at 
times,  the  employment  of  these  standard  values  may  be  reiterated. 
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IV.  Revisiting  the  Results  Presented  in  Reference  1 

It  may  be  useful,  at  this  stage,  to  reproduce  results  that  are  depicted  in 
Reference  1.  To  this  end,  the  following  impositions  are  rendered 

a=>0;  ac=l;  Wcr=0;  (20) 

These  impositions  render  the  complex  commensurate  with  that  defined  in 
Reference  1  and  sketched  in  Fig.  la.  [cf.  Fig.  lb.]  For  these  impositions,  Eqs. 
(16)  and  (18b)  and  Fig.  2a  are  validated.  One  finds  that  Figs.  2a.  1  and  2a.2  are 
akin  to  Figs.  3a  and  7a  of  Reference  1,  respectively.  In  addition  to  evaluating  the 
two-vector  specified  in  Eqs.  (16)  and  (18b),  the  corresponding  induced  reactive 
factor  S^(y)  and  induced  loss  factor  Tj°(y)  are  evaluated  using  Eq.  (8d)  and  the 
assignment  stated  in  Eqs.  (19)  and  (20).  Again,  the  superscript  (o)  in  S(y)  and 
in  7^  (j)  indicates  that  the  evaluations  are  restricted  to  specific  impositions. 
These  evaluations  of  S^(y)  and  of  Tj°(y)  arc  depicted  in  Fig.  3  a.  The  modal 
overlap  parameter  (b)  is  increased  from  the  value  of  (0.1)  to  (2.0)  and  then  onto 
(10)  and  S°(y)  and  Tj^(y)  are  evaluated  and  depicted  in  Figs.  3b  and  3c, 
respectively.  There  are  no  counterparts  to  Figs.  3a.  1,  3b.  1  and  3c.  1  in  Reference 


1.  However,  the  identity  of  Figs.  3a.2  and  3b.2  with  Figs.  5a  and  6a  of  Reference 
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1,  respectively,  is  clear.  In  particular,  the  undulations  that  exist  in  Fig.  3a.2  and 
the  suppression  of  these  undulations  in  Fig.  3b.2  correspond  to  a  phenomenon 
that  is  discussed  in  detail  in  Reference  1.  Figure  3c  does  not  have  a  counterpart  in 
Reference  1.  This  figure  is  included  in  order  to  bring  in  another  phenomenon; 
the  phenomenon  of  erosion  previously  observed  with  respect  to  the  induced  loss 
factor  rjg  (y)  [6,  7].  Again,  the  phenomenon  of  erosion  besets  both,  the  induced 
reactive  factor  iS'(y)  and  the  induced  loss  factor  rj^iy) .  It  transpires  that  mean- 
value  averaging  (arithmatic  mean)  of  the  values  of  S{y)  and  the  mean-value 
averaging  (geometrical  mean)  of  the  levels  of  rj^^y),  respectively,  for  modal 
overlap  parameters  {b)  that  are  small  compared  with  unity,  coincide  and  are  thus 
independent  of  {b)  [1,  17].  This  coincidence  is  illustrated  in  Fig.  3d.  In  this 
figure  S{y)  and  rj^^y)  are  separately  depicted  for  three  small  values  of  the 
modal  overlap  parameter  (Z>);  6  =  (0.01),  (0.1)  and  (0.3).  The  coincident 
curves  representing  the  mean-values  of  S{y)  and  the  mean-levels  of  rj^iy)  for 
all  small  values  of  (Z>);  b  «  (1.0) ,  clearly  emerge  in  Fig.  3d.  These  coincident 
curves  are  designated  5'^(y)  and  rjaiy),  respectively.  An  erosion  is  here  a 
phrase  to  describe  deviations  of  the  values  of  5'(>')  and  of  the  levels  of  7^(7) 
from  these  mean-values  and  mean-levels,  respectively,  when  (b)  approaches  and 
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increases  beyond  the  value  of  unity.  An  erosion  is  then  a  dependence  of  the 
values  of  15(7)  and  of  the  levels  of  77j(y)  on  (b)  as  this  parameter  approaches 
and  exceeds  unity.  To  define  erosion  more  precisely,  a  subsequently  derived  first 
order  approximations  to  .S'(y)  and  to  rigiy)  are  brought  to  bear.  Again,  the  val¬ 
idations  of  the  first  order  approximation  to  the  values  and  levels,  designated 
Sj{y)  and  Tji{y),  respectively,  is  predicated  on  small  values  for  the  modal 
overlap  parameters  {b)\  b«\.  The  curves  that  depict  Sj{y)  and  7/(>’)  are 
dubbed  (FOA).  It  transpires  that  Sj{y)  =  0  and  rij{y)  =  rj^^y) .  Again,  as  the 
values  of  {b)  increase  more  and  more  above  unity,  erosion  in  S{y)  is  manifested 
by  a  progressive  convergence  of  S{y)  from  ■S'^Cy)  onto  5'/(y)  =  0.  On  the  other 
hand,  as  the  values  of  (6)  increase  more  and  more  above  unity,  erosion  in  rj^^y) 
is  manifested  by  a  departure  of  7^(y)fi‘om  7a(y)  =  7/(y)  toward  lower  and 
lower  levels.  Erosions  are  most  pronounced  at  the  edges  of  the  frequency  range 
of  concern;  the  higher  the  values  of  the  modal  overlap  parameters  {b)  above 
unity,  the  more  the  inroad  fi-om  the  edges  into  the  center  of  the  frequency  range. 
(Again,  the  frequency  range  of  reference  is  that  spanning  the  resonance 
frequency  distribution  of  the  satellite  oscillators.)  To  bring  into  focus  the 
existence  and  nature  of  erosion  for  S{y)  and  for  Tj^{y),  as  just  discussed.  Figs. 
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3a,  3b  and  3c  are  overlaid  in  Fig.  3e.  Also,  superimposed  on  Fig.  3e  are  the 
curves  just  dubbed  (FOA).  The  presence  of  erosion,  as  just  described,  is  thus 
revealed  in  Fig.  3e.  This  figure  shows,  however,  that  neither  is  erosion  in  S{y) 
completed  nor  is  erosion  in  Tj^(y)  significant  in  the  vicinity  of  j  =  1 .  To  largely 
complete  the  erosion  in  S(y)  and  to  cause  an  erosion  in  Tjg(y)  at  and  in  the 
vicinity  of  y=l,  the  normalized  overlap  parameter  (b);  b  =  [(;r6)(i?+l)“*  ], 
needs  approach  and  exceeds  unity  [1,  6,  7].  This  extreme  case  of  erosion  is 
illustrated  in  Fig.  3f.  In  this  figure  the  only  change,  in  parameters  that  specify 
the  complex  assigned  to  Fig.  3e,  is  the  number  (R)  of  satellite  oscillators.  The 
number  (R)  is  (7)  instead  of  (27),  so  that  for  b  ==  10,  (6)  comfortably  exceeds 
unity  in  Fig.  3f  After  this  brief  digression,  it  is  time  to  return  to  consideration  of 
Reference  1  and  beyond. 

The  limited  scope  of  Reference  1  curtails  the  modeling  and  the  analysis  of  the 
complex  there  considered,  notwithstanding  that  complexes  employed  to  date  are 
largely  subjected  to  similar  limitations  [2-1 1].  A  poignant  question  arises:  what 
is  fi'esh  about  the  complex  defined  herein  as  compared  with  the  complex  defined 
in  Reference  1?  Whereas  in  Reference  1  a  satellite  oscillator  in  isolation  is 


characterized  by  a  mere  mass  control  term,  here  it  is  characterized  by  the 
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oscillator  impedance  Z^(ft)),  as  stated  in  Eq.  (3b).  This  oscillator,  in  addition  to 
the  mass  control  term,  may  also  possess  a  stiffiiess  control  term.  Moreover,  the 
coupling  impedance  as  stated  in  Eq.  (3c),  may,  in  addition  to  the 

stiffiiess  control  term,  possess  also  a  mass  control  term.  Finally,  the  coupling 
between  a  satellite  oscillator  and  a  master  oscillator  may  be  allowed  to  include  a 
gyroscopic  control  term  [12,  13].  Obviously,  the  complex  sketched  in  Fig.  lb 
and  formulated  in  Eqs.  (5)  -  (8)  is  more  versatile  than  that  in  Reference  1,  as 
sketched  in  Fig.  la.  It  may  be  useful,  therefore,  to  investigate  a  few  of  the 
attributes  of  this  versatility  even  under  the  similarity  conditions  imposed  in  Eq. 
(11)  and  the  simplifying  assumptions  proposed  in  Eqs.  (14)  and  (15). 


V.  Various  Coupling  Forms  and  Coupling  Strengths 


It  is  of  interest  to  compute  and  display  the  induced  reactive  factor  S(y)  and 
the  induced  loss  factor  ,  as  functions  of  (y),  for  a  variety  of  selected  coup¬ 
ling  forms  and  coupling  strengths.  The  coupling  forms  are  defined  according  to 
whether  the  coupling  is  dominated  by  either  a  stiffiiess  coupling,  a  gyroscopic 
coupling,  a  mass  coupling  or  combinations  of  these  coupling  types.  The 
coupling  strengths  of  these  forms  are  determined  by  the  values  of  the  coupling 
parameters;  the  stiffness  coupling  parameter  (a^),  the  gyroscopic  coupling 
parameter  (g);  gr  =g,  and  the  mass  coupling  parameter  (m^) .  The  values  of 
these  coupling  parameters  may  be  categorized  from  weaker-coupling,  to 
moderate-coupling  onto  stronger-coupling  in  the  range  of  values  of  0.03,  0.15 
and  0.75,  respectively.  In  this  categorization,  the  coupling  of  the  satellite 
oscillators  to  the  master  oscillator,  defined  in  Reference  1,  is  of  stiffiiess  control 
coupling  form;  i.e.,  ^  0,  g  =  mc  =  0,  and  of  (very)  strong  coupling  strength; 

namely,  =  1.0  [a =0.0.].  Such  a  coupling  form  and  a  coupling  strength  define 
satellite  oscillators  commonly  designated  sprung-masses  [1-11].  The  new 
evaluations  in  this  paper  are  exhibited  in  Figs.  4-7.  Of  significance  are  not  only 
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the  variations  in  the  coupling  forms  and  in  the  coupling  strengths,  but,  also,  the 
influence  that  changes  in  the  modal  overlap  parameters  have  on  the  nature  of  the 
induced  reactive  factor  5'(y)  and  induced  loss  factor  77^(7) .  A  first  set  of  figures 
is  evaluated  with  ^>  =  (0.1) ,  a  second  with  Z>=(2.0)  and  a  third  with  6=(10) .  A 
major  feature,  common  to  all  evaluations,  is  that  the  undulations  in  the  first  set, 
for  which  6 =(0.1),  is  suppressed  in  the  second  and  in  the  third,  for  which 
6=(2.0)  and  6=(10),  respectively  [1].  Figures  4-7  are  cast  in  the  format  of  Fig. 
3e;  corresponding  figures  in  the  three  sets  are  overlaid  so  that  the  undulations  in 
the  first  set  and  their  suppression  in  the  second  and  third  are  observed  at  a  glance. 
Another  major  feature,  common  to  all  evaluations,  is  that  the  erosions,  discussed 
briefly  in  the  preceding  Section  with  respect  to  S(y)  and  7^(3;)  and  depicted  in 
Figs.  3e  and  3f,  thread  all  figures.  It  is  noted  that  when  (b)  [=  (7rb)(R+l)~^  ] 
approaches  and  exceeds  unity,  erosions  in  S(y)  are  nearly  completed  and  in 
77^(3;)  erosions  occur  even  at  and  in  the  vicinity  of  3^=1.  These  extreme 
erosions,  however,  in  addition  to  the  dependence  on  (6),  seem  to  carry  a  slight 
dependence  on  the  coupling  forms  and  strengths.  These  dependencies,  which  are 
more  clearly  apparent  when  the  coupling  is  weak;  e.g.,  in  Figs.  6b  and  7d,  are,  at 


this  stage,  merely  noted,  notwithstanding  that  the  linear  scale,  reserved  for  the 
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induced  reactive  factor  S(y),  hardly  exhibits  values  that  belong  to  moderate  and 
weak  coupling  strengths.  Supplementally,  some  of  the  emerging  details  in  Figs. 
4-7  may  be  summarized  as  follows: 

1.  The  format  of  Fig.  3e  is  reproduced  in  Fig.  4a,  except  that  the  stiffiiess 
coupling  parameter  is  reduced  from  =  1.0  to  =  0.75 ;  i.e.,  from  very  strong 
to  strong  stiffiiess  control  coupling  .  This  change  decreases  the  values  of  S(y) 
and  the  levels  of  7^(y)in  Fig.  4a  as  compared  with  those  in  Fig.  3e.  The 
decrease  is,  however,  slight.  A  more  drastic  decrease  in  the  values  of  S(y)  and 
in  the  levels  of  rjg(y)  occurs  in  Fig.  4b  as  compared  with  Fig.  3e.  In  Fig.  4b  the 
stiffiiess  coupling  parameter  is  a^=0.15,  whereas  in  Fig.  3e  «^=1.0.  This 
decrease  in  the  values  of  S{y)  and  in  the  levels  of  J]^(y)  is  largely  related  to  the 
difference  in  the  coupling  strengths.  In  Fig.  4b  the  coupling  strength  is  moderate. 
Again,  it  is  noted  that  on  the  linear  scale,  on  which  the  induced  reactive  factor 
S(y)  is  a  priori  depicted,  the  values,  even  for  moderate  coupling  strengths,  are 
well  nigh  negligible.  (There  is  little  logic  to  expanding  the  linear  scale  to  gain 
insight  into  the  behavior  of  negligible  values.  The  negligibility  itself  is  infor¬ 


mation  enough.) 
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2.  The  format  of  Figs.  4a  and  4b  is  reproduced  in  Figs.  5a  and  5b,  respectively, 
except  that  the  stiffiiess  coupling  form  is  changed  to  a  gyroscopic  coupling  form; 
namely,  =  m^=0,  g^O.  In  Fig.  5a  the  coupling  is  strong;  g=0.75 ,  and  in 
Fig.  5b  the  coupling  is  moderate;  g=0.15 .  The  similarity  between  Figs.  4  and  5 
is  obvious.  Also  obvious  is  the  slope  in  the  curves  in  Fig.  5.  This  slope  is 
characteristic  of  the  gyroscopic  coupling.  The  gyroscopic  coupling  enters  in  the 
form  of  (g/y)  and  not  merely  in  the  form  of  (g) .  [cf  Eq.  (8),] 

3.  The  format  of  Figs.  4a  and  4b  is  reproduced  in  Figs.  6a  and  6b,  respectively, 

except  that  the  stifftiess  coupling  form  is  changed  to  a  mass  coupling  form; 
namely,  a^=g=0,  0.  In  Fig.  6a  the  coupling  is  moderate;  w^=0.15,  and 

in  Fig.  6b  the  coupling  is  weak;  m^,  =0.03 .  The  similarity  between  Figs.  4  and  6 
is  obvious  although  the  values  of  .^(y)  and  the  levels  of  in  the  former  are 
higher  than  in  the  latter,  largely  in  consequence  of  the  disparities  in  the  coupling 
strengths. 

4.  Finally,  the  format  of  Fig.  4a  is  reproduced  in  Fig.  7,  except  that  the  stiffness 
coupling  form  is  modified  to  accommodate,  in  addition,  another  form  of 
coupling.  In  Fig.  7a,  the  additional  coupling  form  is  mass  control;  namely: 

=  Wc=0.75[a=0.1.],  g=0,  which  is  a  very  strong  coupling  strength.  In 
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Figs.  7b,  7c  and  7d,  the  additional  coupling  is  gyroscopic  control;  namely,  in  Fig. 
7b:  0.53  [a =0.43.],  g=0.54,  fn^  =  0,  which  is  a  strong  coupling  strength, 

in  Fig.  7c:  =  0.10  [a =0.9.],  g  =  0.1 1,  =  0,  which  is  a  moderate  coupling 

strength  and  finally,  in  Fig.  7d:  =  0.02  [ar=0.98.],  g  =  0.022,  =0,  which 

is  commensurate  with  a  weak  coupling  strength.  The  values  of  S{y)  and  the 
levels  of  rjgiy)  in  these  figures  are  set  largely  by  the  coupling  strengths;  e.g.,  the 
higher  the  coupling  strength  the  higher  are  the  levels  of  ri^{y). 

In  the  normalized  frequency  range  of  concern,  the  values  of  S{y)  and  the  levels 
of  ri^{y)  in  Figs.  3-7  that  pertain  to  modal  overlap  parameters  (b)  that  are  small 
compared  with  unity;  e.g.,  b  =  (0.1),  undulate.  On  the  other  hand,  again  in  the 
normalized  frequency  range  of  concern,  the  curves  in  Figs.  3-7  that  pertain  to 
modal  overlap  parameters  that  exceed  unity;  e.g.,  Z>=(2.0)and  (10),  are 
reasonably  smooth.  A  few  questions  arise:  Can  these  features  be  estimated  by 
replacing  the  summation  in  Eq.  (8d)  by  an  integration  and  if  so,  can  this  integral 
be  performed  with  ease?  Can  the  result  of  this  performance  interpret  the 
response  behavior  of  the  master  oscillator  in  terms  of  the  various  parameters  that 
define  the  complex?  And  last,  but  by  no  means  least,  what  about  the  undu¬ 


lations,  do  they  feature  in  the  result  of  this  integration? 
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VI.  Replacing  a  Summation  by  an  Integration 

The  index  (r)  is  given  a  continuous  connotation  and  the  summation  in  Eq. 

(8b)  is  replaced  by  integration.  Under  a  condition  that  allows  this  replacement, 

Eq.  (8b)  assumes  the  form  [1] 

,  rz^(M+s) 

[S (y)  -  i  ri^  (y)]  =  (y)^  j  dz^  {r  )  [/(/*)  //(/*)]• 

\(,e) 

{[1  -  {z{r)f  [fncir)-{Zc(f)f  (1  +  //7c(^))]  “  }  • 

{[l+fh^(f)\-[{Zr{r)f  (l+/;7r(^))]^^  (21) 

where 

mr)lyf  =  4m,(F)  {z,{r)f  {Uirj.if))  +Kg{r)lyf\ 

m{r)= /i{r){R  +  \y^ f(7)dz^(7)  =  dF;  {z^(r)}^=  {z(F)}^+  {z^(F)}^; 

F  =  f(/?+l)"*;  (^+F)  =  (/?+f)(i?+l)"^  (22) 

If  is  small  enough  so  that  the  vanishing  of  the  real  part  of  the  denominator 
in  the  integral  in  Eq.  (20)  predominates  the  values  of  the  integral,  the  integral 
yields 

)}^  7(j^)  -  {^C  (^)}^  7c  (^  )]+ 4  (F)  (z^  (7^  )}^  7c  (7^ )} 
-i[mc(ro)-{ZcO"o)ff  +  [M^)^c(^)}^7(^)7c(^)]+[7c(^)/Ff }].  (23) 
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where 

(J )  <  >^  <  (^ ); 

=(l  +  ^cr)  =  {2(^)}^  +  {Zc{ro)f', 

{\  (^)}^  7r„ (^)={2(^)}^  V{ro)  +  {Zciro)f  Vd^o)  •  (24) 

Equation  (24)  defines  the  specific  value  of  (^)  and  it  is  recognized  that  (^)  is  a 
function  of  (y)  and  vice  versa.  Adopting  the  impositions  and  the  simplifications 
that  are  conducted  in  Eqs.  (16)  and  (18),  one  may  derive  fi*om  Eqs.  (23)  and  (24) 
the  results 

S{y)  =  DA{n{y)y,  S{y)=^SM  =  0,  (25a) 

^siy)=I^{C  +  0{n{y)f-y,  Ti^{y)^7]j{y)  =  DC,  (25b) 

with 

D  =  [TuHyilmMs  mK)]  [(l+m,)]-' ,  (26a) 

A  =  {{fhc  +  a^f  +  {nic  -  ayi\,  (26b) 

C  =  [{fhc  +ac)^  +iglyf] ,  (26c) 

O  =  [(1+ ] ,  (26d) 

and  where 

/('i)=((a+«c/''^V(«)W’r‘;  /<(ii)=(Jl/,/Mfl),  (26e) 

';(:>')  = '7('i)  =  (f/®)I>W(y)^];  fif/i)  =  tA'i)  =  \<ro)>  (26f) 
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(1+m.)  =  (a+a,) ;  [l+(?'/2)r'‘'^>  <  3-  £  [l-(y/2)r<''^> .  (26g) 

The  underlying  condition  of  the  validity  on  Eq.  (25)  is  detailed  in  Reference  1. 
This  condition  holds  even  though  the  coupling  forms  are  elaborated  to  include 
not  only  a  weaker  stiffhess  coupling,  but  also  mass  and  gyroscopic  coupling 
forms.  Moreover,  the  coupling  parameters  may  define  various  degrees  of 
coupling  strengths.  Strictly,  the  validity  of  Eq.  (25)  demands  the  equality  of 
S(y)  to  the  primary  term  iS'/(32)and  75(y)to  the  primary  term  T]j{y).  These 
equalities  designate,  ^/(y)  and  rjjiy)  the  first  order  approximations  to iS'(>’)  and 
to  ^^(y),  respectively.  The  terms  A{T]{y)}  and  ^{^(j^)}^  are  of  the  order  of  the 
higher  approximations  to  the  integral,  notwithstanding  that  situations  exist  in 
which  (A)  and  (O)  are  identically  equal  to  zero;  e.g.,  when  the  satellite 
oscillators  are  sprung-masses  for  which  =1,  g"  =  =  0.  [cf  Eq.  (20).]  In 

these  situations  the  equality  of  the  integral  evaluations  of  S(y)  to  S/(y)  and  of 
fjjiy)  need  not  be  specifically  invoked.  (More  on  this  subject  when  in 
a  subsequent  paper  higher  order  approximations,  than  the  first,  are  to  be 
evaluated.)  Clearly  and  significantly,  the  primary  terms  Sf(y)  and  p/(y)  and, 
therefore,  the  first  order  approximations  of  iS'(>')  and  of  %(y)  are  independent 


of  the  modal  overlap  parameter  (b).  Without  much-a-do,  in  this  paper  only  the 
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first  order  approximations  are  implemented  and  considered.  In  this  approx¬ 
imation  the  first  order  approximation  to  S{y)  is  the  primary  term  Sj{y)  and  to 
risiy)  is  the  primary  term  tJiiy) . 

Equation  (25)  reveals  the  parametric  composition  of  the  first  order  approx¬ 
imation  to  the  induced  reactive  factor  5'/(y)  and  to  the  induced  loss 

factor  rig{y)  T]j{y).  Again,  one  is  reminded  that  the  induced  reactive  factor 
S{y)  and  the  induced  loss  factor  describe  the  influence  of  the  coupled 
satellite  oscillators  on  the  impedance  of  the  master  oscillator.  In  the  absence  of 
coupling  and  in  units  of  the  reactive  factor  is  [j^-1];  in  the 

presence  of  coupling  it  is  [y^  +  .  On  the  other  hand,  in  the  absence  of 

couplings  the  loss  factor  is  (7^);  in  the  presence  of  couplings  it  is  [7o+7j  (>’)]. 
When  S{y)  is  positive  the  satellite  oscillators  contribute,  through  the  couplings, 
a  mass-like  term,  when  negative,  a  stiffiiess-like  term.  On  the  other  hand,  since 
the  satellite  oscillators  add  merely  passive  elements  to  the  complex,  7^(7)  is 
invariably  positive.  Equation  (25)  confirms  this  statement  and  exposes  the 
proportionality  of  Sj{y)  and  of  7/(7)to  (D).  Therefore,  Sj{y)  and  rijiy)  are 
directly  proportional  to  the  mass  ratio  (M_j/M);  it  is  estimated  that  for  a 
reasonable  complex  with  !  M)  equal  to  about  a  tenth,  (D)  is  of  the  order  of 
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unity,  [cf.  Eq,  (19).]  The  equality  of  Sj(y)  to  zero  and  the  quadratic 
dependence  of  the  primary  term  Tjjiy)  of  Tjgiy),  in  terms  of  (C),  which 
entertain  the  term-components  and  (g/y)^,  is  of  significance. 

Again,  it  is  emphasized  that  the  primary  terms  Sj(y)  and  Tjj(y)arQ  the  true  first 
order  approximation  to  S(y)  and  to  ?Js{y),  respectively;  in  this  context  the 
terms  DA  {rjiy)}  and  DO{ri{y)f'  are  superfluous.  In  any  case,  Eq.  (26) 
indicates  that  DO  {niy))  in  rj^iy)  is  rarely  dominant  even  when  the  loss  factor 
T]{y)  of  a  typical  satellite  oscillator  exceeds  any  of  the  coupling  parameters; 

y)'^'niy)y  notwithstanding  that  in  the  absence  of  any  couplings 
^'(y)  and  rj^{y) ,  as  stated  in  Eq.  (25),  are  negligible  on  account  of  .<4=0  and  of 
C=0  and  0  =  0,  respectively.  When  couplings  do  exist,  there  are  values  of 
S{y)  that  are  non-zero  whereas  those  of  Sj{y)  are  equal  to  zero.  Yet  when  the 
coupling  strengths  are  weak,  and  even  moderate,  the  values  of  |  <5(7)  |  are  small 
compared  with  unity  even  if  (y)  is  at  the  edges  of  the  frequency  range  of 
concern.  On  the  other  hand,  when  couplings  do  exist,  the  levels  in  both,  the  first 
order  approximation  and  the  exact  evaluations  of  77j(y),  are  certainly  not 
negligible,  even  if  the  couplings  are  weak.  (Note  that  (7/0) of  the  order  of 
(10“^)is  not  unreasonable  [12].) 
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As  already  discussed  and  demonstrated  in  Figs.  3-7  the  curves  for  a  modal 
overlap  parameter  (b)  that  is  small  compared  with  unity;  b«l,  possess  values 
and  levels  that  are  undulated.  The  excursions  in  the  undulations  are  the  more 
pronounced  the  smaller  is  the  value  of  (b)  [11].  Nonetheless,  the  mean-value 
averaging  of  the  undulating  values  of  S{y)  and  of  the  vmdulating  levels  of  Tjg{y) 
converge  unto  two  single  curves,  which  are  designated  S^iy)  and  ?Ja(y)’  respec¬ 
tively.  [cf.  Fig.  3d.]  As  (Z>)  increases,  approaching  and  exceeding  unity,  the 
undulations,  again  true  to  form,  are  suppressed.  As  already  intimated  in  Section 
III,  increasing  (6)  beyond  unity  brings  in  the  phenomenon  of  erosion  which 
increases  as  the  modal  overlap  parameter  (b)  reaches  higher  and  higher  above 
unity. 

To  accentuate  the  nature  of  erosion  in  S(y)  and  in  Tig(y),  the  first  order 
approximation  (FOA)  of  S(y)  [=5/(7)]  and  of  ^j(7)[=77/(7)]  are  appropriately 
superimposed  on  Figs.  4-7.  [cf.  Figs.  3d,  3e  and  3f.]  It  is  now  observed,  in  these 
figures,  that  the  mean-value  averaging  of  the  undulations  of  the  exact  values  of 
S(y)  and  of  the  exact  levels  of  rj^iy) ,  when  (b)  is  small  compared  with  unity; 
b«l,  relate  to  the  first  order  approximation  (FOA)  of  5(7)  and  of  7^(7)  in 
the  form:  Sa(y)^Sj{y)  except  at  7=1  and  7^(7)  =  7/(7),  respectively.  Since 


46 


Sa(y)  and  Sj(y)  and  Tj^iy)  and  7jj(y)  are  independent  of  (b),  these  relation¬ 
ships  are  also  independent  of  (b).  [cf  Figs.  3d  and  3e.]  (The  independence  of 
7/(t)  of  (b)  has  been  stretched  by  some  to  conclude  that(Z))  may  be  rendered,  a 
priori,  equal  to  zero.  Neglecting  to  mention  in  this  rendering  that  mean-value 
averaged  levels  are  substituted  for  highly  undulated  levels,  is  not  a  viable 
scientific  procedure,  unless  ignorance  is  bliss  [11].)  On  the  other  hand,  when  (b) 
approaches  and  exceeds  unity,  the  exact  values  of  S(y)  and  the  exact  levels  of 
f/s(y)  become  free  of  undulations,  but  these  values  and  levels  erode  with  further 
increases  of  (b)  [1,  6,  7,  15].  The  erosion  in  the  exact  values  of  <5(7)  is  a 
departure  from  the  curve  designated  .5^(7),  progressively  converging  on  the  first 
order  approximation  at  Sj(y)  =  0,  as  (b)  increases  higher  and  higher  above 
unity.  In  contrast,  the  erosion  in  the  exact  levels  of  commences  and 

progresses  from  the  levels  of  the  first  order  approximation,  at  7/(y),  toward 
lower  and  lower  levels,  as  (b)  increases  higher  and  higher  above  unity.  To 
account  for  these  progressive  dependencies  on  (b) ,  higher  order  approximations 
are  clearly  required  [20].  What  is  doubtful  is  whether  higher  and  higher  order 
approximations  can  account  for  the  undulations  when  (b)  needs  to  remain  small 
compared  with  unity.  To  account  for  the  undulations,  an  entirely  different 
approximation  procedure  is  thus  called  for.  In  this  paper  devising  such  an 
approximation  procedure  is  not  attempted. 


VII.  A  Typical  Member  of  an  Ensemble  of  Complexes  Supporting 
Various  Parametric  Combinations 

In  the  preceding  evaluations,  the  distribution  of  resonance  frequencies 
and  the  assigned  loss  factors  (77^^)  for  the  satellite  oscillators  are  sequen¬ 
tial  functions  of  the  normalized  index  (F) .  These  two  quantities,  exemplified  in 
Figs.  2a.l,  2b.l  and  2c.l  and  in  Figs.  2a.2,  2b.2  and  2c.2,  respectively,  may  be 
smoothed  out  by  extrapolation  and  interpolation  into  monotonic  and  continuous 
ftinctions  of  (F);  (F)=[r(J?  +  l)“^  ].  This  kind  of  smoothness  is  rarely  found  in 
practice  and  a  question  arises  as  to  what  are  the  expected  consequences  of  more 
practical  assignments  for  these  parameters  and  others?  In  this  section  a  few 
layers  are  removed  in  the  quest  to  discover  the  phenomena  that  may  be 
encountered,  in  the  induced  reactive  factor  iS'(j^)  and  in  the  induced  loss  factor 
rig{y),  by  the  insertion  of  these  more  realistic  parametric  values.  Since  the 
assignment  for  die  parameters  that  define  the  satellite  oscillators  and  their 
couplings  to  the  master  oscillator  can  hardly  be  drawn,  a  more  generalized 
approach  is  undertaken  to  investigate  the  influence  of  introducing  variations  in 
these  parametric  values.  In  particular,  in  this  section  two  parameters  are  selected 
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to  carry  these  variations;  either  individually  or  in  unison.  In  the  first,  the  index 
(r)  of  a  satellite  oscillator  is  assigned  a  distribution  of  pseudo-statistical  values. 
[Pseudo-statistical  is  in  reference  to  a  sample  selected  out  of  an  ensemble  of 
random  samples.]  The  index  (r)  is  distributed  sequentially  and  fractionally,  in 
the  range  1  <  r  <  27.  A  pseudo-statistical  index  is  designated  A(r) ,  where 
A(a*)  <  A(q);  q  =  (r  +  l);  \<r  <(R-l).  [cf.  Eq.  (9b).]  In  the  second,  the 
modal  overlap  parameter  (6;.)  is  assigned  a  distribution  of  pseudo-statistical 
values  that  span  the  ranges  (2)  >  >  (0.1)  and  (3.5)  >b^>  (1),  respectively. 

The  distribution  of  A(r)  and  of  (b^),  with i?  =  27,  that  are  employed  in  this 
section  are  depicted  graphically  in  Figs.  8a,  8b,  and  8c.  The  two-vector 
{(^/•r)»  (7rr)}(^)»  Stated  in  Eqs.  (16)  and  (18),  is  typically  depicted,  for  the 
pseudo-statistical  values  shown  in  Figs.  8a,  8b,  and  8c,  in  Figs.  9a,  9b,  and  9c, 
respectively.  Figure  9a  depicts  (jc^^)  as  a  function  of  A{r)  and  Figs.  9b  and  9c 
depict  (fJrr)’  ^  ^  function  of  (r),  where  A(r)  =  [A(r)(7?+1)“*]  and 
r  =  [r(ii!  +  l)“*].  [cf.  Fig.  2a.]  It  is  observed,  in  Fig.  9a,  that  the  pseudo- 
statistical  variations  embody  the  phenomenon  of  mode  bunching  in  which 
variations  in  the  modal  density  of  the  satellite  oscillators  drastically  vary  as  a 
function  of  A(r)  [21].  On  the  other  hand,  as  Figs.  9b  and  9c  show,  the  loss 
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factor  as  a  function  of  (r),  faithfully  follows  the  variations  assigned  to 

ipr) .  In  Fig.  9b  some  of  the  values  of  (Z>^)  are  less  than  unity,  in  Fig.  9c  all  the 
values  of  (Z>^)  are  in  excess  of  unity. 

The  influence  of  the  variations,  described  in  Fig.  8,  on  the  induced 
reactive  factor  S{y)  and  on  the  induced  loss  factor  %(y),  as  functions  of  (y), 
are  exemplified  in  Figs.  10.1  10.2,  as  a  pair,  and  in  Figs.  11.1  and  1 1.2.,  as  a  pair, 
respectively.  Each  pair  of  figures  represents  a  set  of  figures.  The  first  pair  of 
figures  in  each  set;  e.g..  Figs.  10.1a  and  10.2a  depict  the  base  situation  in  which 
A(r)  =  F  and  =  \.  The  second  pair  of  figures  in  each  set;  e.g..  Figs.  10.1b 
and  10.2b,  depict  the  situation  in  which  A(r)  is  as  shown  in  Fig.  8a  and  =1 . 
The  third  pair  of  figures  in  each  set;  e.g..  Fig.  10.1c  and  10.2c,  depicts  the 
situation  in  which  A(r)  =  r  and  (6^)  is  as  shown  in  Fig.  8b.  The  fourth  pair  of 
figures  in  each  set;  e.g..  Figs  10.  Id  and  10.2d,  depict  the  situation  in  which 
A(r)  =  F  and  (b^)  is  as  shown  in  Fig.  8c.  The  fifth  pair  of  figures  in  each  set; 
e.g..  Figs.  10. le  and  10.2e,  depict  the  combined  situation  in  which  A(r)  and 
ibj.)  are  as  shown  in  Figs.  8a  and  8b,  respectively.  Finally,  the  sixth  pair  of 
figures  in  each  set;  e.g..  Fig.  10. If  and  10.2f,  depict  the  combined  situation  in 
which  A(r)  and  (Z>^)  are  as  shown  in  Figs.  8a  and  8c,  respectively.  It  is 
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recognized  then  that  each  pair  of  figures  in  a  set  presents  a  complete  evolution  in 
the  process  of  applying  the  pseudo-statistical  variations  depicted  in  Fig.  8  to  the 
two  parameters  A(r)  and  {b^).  Also,  each  pair  of  figures  in  a  set  selects  a 
specific  coupling  form  and  a  specific  coupling  strength.  Thus,  Fig.  10  pertains  to 
a  strong  stiffiiess  coupling:  =1.0  [a =0.0.],  =  f  =  0,  and  Fig.  1 1  pertains 

to  a  mix  of  stiffiiess  and  of  gyroscopic  coupling  of  moderate  strength: 

=  0.1  [«  =  O.9.],  g  =  0.1 1,  =  0.  [cf  Figs.  3e  and  7c.]  [One  recalls  that  the 

values  of  S{y)  relating  to  moderate  and  especially  weak  coupling  strengths  are 
largely  negligible  compared  with  unity,  even  if  (b)  is  as  small  as  (0.1). 
Therefore,  the  values  of  S(y)  depicted  in  Fig.  1 1  are  hardly  discernible.] 

The  first  pair  of  figures  of  each  set;  namely.  Figs.  10.1a  and  10.2a,  and  1 1.1a 
and  1 1 .2a  exhibit  undulations  in  the  values  of  the  induced  reactive  factor  S(y) 
and  in  the  levels  of  the  induced  loss  factor  Tjg(y) ,  as  functions  of  (y).  However, 
these  undulations  are  small  and  they  are  completely  suppressed  as  soon  as  (b^) 
approaches  the  value  of  (2).  [cf  Fig.  3b.]  The  variations  depicted  in  Fig.  8  are 
clearly  discernible  in  all  the  subsequent  figures  in  the  series  entitled  Figs.  10  and 
11.  True  to  form,  there  is  but  a  tinge  of  edge  erosion  in  Figs.  10.1a  and  10.2a, 


and  in  Figs.  11.1a  and  11.2a  yet  in  the  likes  of  Fig.  3b  stronger  sign  of  erosion 
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has  already  reared  its  head.  From  Figs  10.1a  and  10.2a,  and  jfrom  Figs.  11.1a  and 
11.2a  to  Fig.  3b,  (b)  is  changed  from  unity  to  merely  two.  To  confirm  this 
statement  and  to  provide  for  convenient  and  interpretable  data  from  which  to 
judge  the  more  erratic  data  that  incorporates  the  pseudo-statistical  variations,  the 
first  order  approximation  of  S{y)  and  of  ?]s(y),  given  in  Eq.  (25),  are  super¬ 
imposed  on  Figs.  10.1a  and  10.2a,  and  on  Figs.  11.1a  and  11.2a  and  on  all  other 
figures  in  the  series  entitled  Figs.  10  and  11. 

Again,  the  pseudo-statistical  variations  are  defined  by  two  competing  and 
nearly  independent  factors,  i.e.,  by  A(r)  and  ^(a*)[=(Z>^)].  In  Figs.  9a,  10.1b, 
10.2b,  11.1b  and  11.2b;  it  is  observed  that  at  a  mode  bunching  (a  rich  modal 
density)  region  the  influence  of  the  satellite  oscillators  is  more  pronounced  than 
at  a  mode  sparsity  (a  poor  modal  density)  region  [21].  On  the  other  hand,  when 
the  modal  overlap  parameter  (6^)  entertains  values  that  are  small  compared  with 
unity,  the  values  and  levels,  as  functions  of  (y),  tend  to  fluctuate.  The  fluct¬ 
uations  are  pronounced  at  and  in  the  vicinity  of  the  resonance  frequencies  of 
those  satellite  oscillators  to  which  these  small  values  of  (6^)  are  assigned.  At  and 
in  the  vicinity  of  the  resonance  frequencies  of  those  satellite  oscillators  to  which 


(b^)  are  assigned  values  that  approach  and  exceed  unity,  no  such  fluctuations  are 
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present;  e.g.,  see  Figs.  10.1c,  10.2c,  11.1c  and  11.2c,  and  Figs.  lO.ld,  10.2c, 
11. Id  and  11. 2d  and  contrast  them,  respectively  [22].  When  variations  in  both 
parameters  are  combined,  both  characteristics  can  be  identified  in  the  values  of 
the  induced  reactive  factor  S(y)  and  in  the  levels  of  the  induced  loss  factor 
Vs(y)’  see  Figs.  10  and  11  and  contrast,  in  particular,  Figs.  lOle,  10.2e, 
1  Lie  and  1 1.2e  with  Figs.  10. If,  10.2f,  11. If  and  11.2f,  respectively. 
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Appendix  A 

The  sweeping  assumptions  rendered  in  Eq.  (15),  which  leads  to  Eqs.  (9d) 
and  (16),  culminating  in  Eq.  (18),  may  be  introduced  more  gradually.  The 
purpose  for  this  Appendix  is  to  effect  such  a  gradual  introduction.  In  this  maimer 
when  some  of  these  assumptions  are  relieved,  reevaluation  of  the  induced 
reactive  factor  S(y)  and  of  the  induced  loss  factor  ?Js(y)  may  be  readily  insti¬ 
tuted. 

From  Eqs.  (9c)  and  (12),  the  design  demands  that 
(«r  +  «cr)  for  F=(l/2),  (Al) 

and  if  (x°)  becomes,  by  design,  unity  at  r  =(1/2),  the  expression  reduces  to 
(ttr  +acr)  =  (1+ ^cr)  ^r  F  =  (1/2).  (A2a) 

From  Eq.  (12)  it  is  observed  with  satisfaction  that  (x°)  is,  indeed,  unity  at 
F  =  (1/2) .  On  the  other  hand,  from  Eqs.  (10)  and  (11)  one  obtains 
(x,fnr=<.brl2)[d{x,(,x;f}/er];  (x^fticr  =(*cr/2)[3{a„(*?)^}/ar]; 

=[2(a.+a„)(x,“)^r'  {Md{arix°f}ldr]*b„[d{a„(x°f)ldr-\},  (A3) 
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where  again,  (r)  is  allowed  to  have  a  continuous  connotation  as  explained  in 
Reference  1.  In  particular,  if  the  spring  factors  (a^)  and  (a^r)  independent 
of  (r) ,  Eq.  (A3)  simplifies  to  read 
>7,  =  (b, ];  =  (b,,  /^)[r(R)(x^f  ]; 

nrr  =  PrVcr'^  CL\  a^r  =  »  (A4) 

where  (x^)  is  extrapolated  and  interpolated  to  become  a  continuous  function 
of(r) 

[Oln(x»/ar)]=  («  +  !)-■  [/(RXx^)^],  (A5) 

(A-  A;r)  =  (a  +  «£>“'  ((*r  “c)].  (A6) 

b,=[(;rb,)(Jt  +  l)-'];  A„=[(;rft„)(S  +  l)-'],  (A7) 

and  from  Eq.  (A2a),  by  design 

(a+(Zc)  =  (l+m^).  (A2b) 

From  Eqs.  (10)-(12)  and  (A2)-(A7)  one  may  cast  the  designed  two- vector  in  the 
form 

{(x„),  (>;,,)}(«)  =  +  .  {(x»  ),  (>?«))  (R),  (A8) 

(K),  «))  W  =  {(*?).  />^)lr(R)(x^f  ]},  (A9) 


where 
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i\+mc)'^[a{bj.lbcr)+ac\  ;  b^^b^,,  (AlOa) 

(l  +  mc)”^[a(^/6c)+aj;  b^=b,  bcr  =  bc  (AlOb) 

=  P=< 

1  ;  b  =  br  =  bcr,  (AlOc) 

the  quantities  (m^),  (x^),  y(R)  and  (b^^)  are  stated  in  Eqs.  (Id),  (12),  (13)  and 
(A7),  respectively,  and  the  superscript  (o)  in  (x^^)  and  (7^^)  recognizes  that 
(x^^)  and  (7®^)  are  restricted  to  specific  impositions,  [cf.  Eq.  (4).]  For  the  sake 
of  computational  and  interpretive  advantage  the  validity  of  Eqs.  (A4)  and  (A  10c) 
is  universally  adopted  in  this  paper,  [cf.  Eq.  (15).]  Under  this  imposition,  Eqs. 
(A8)  and  (A4)  simplify 

=  +  A=[(^A)(J?  +  ir' ].  (All) 

tjCr)  =  rtr=  Vcr  =  Vrr  S  '7('^)  =  (*  ! ’t)[r(R)(x°f\  (A12) 

where  (x°)  is  stated  in  Eq.  (12),  and  (J3^)  becomes  equal  to  unity,  [cf.  Eqs. 


(18a)  and  (16).] 
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Fig.  la.  A  master  oscillator  attached  to  a  set  of  sprung-masses. 


Fig.  lb.  A  set  of  satellite  oscillators  coupled  to  a  master  oscillator. 


Fig.  2.  (1.)  Resonance  frequency  distribution  and  (2.)  Corresponding  as¬ 
signed  loss  factor  (;7;.^)for  the  satellite  oscillators  as  a  function  of  the  normalized 
index  (F) ,  respectively.  (*) 

a.l  and  a.2.  With  R  =  27  and  =  0.  ,  t],,  =  ri% .] 

(*)  When  regions  of  curves  clearly  overlap,  the  color  of  the  one  with  the  higher  modal 
overlap  parameter  (b)  wins. 
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Fig.  2.  (1.)  Resonance  frequency  distribution  (x^^)  and  (2.)  Corresponding  as¬ 
signed  loss  factor  (77^^)  for  the  satellite  oscillators  as  a  function  of  the  normalized 


index  (r),  respectively.  (*) 

b.l  and  b.2.  As  in  a.l  and  a.2,  exeept  that  (R)  is  changed  from  (27)  to  (7). 
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Fig.  2.  (1.)  Resonance  frequency  distribution  (x,,)  and  (2.)  Corresponding  as¬ 
signed  loss  factor  (;7^;.)for  the  satellite  oscillators  as  a  function  of  the  normalized 
index  (r),  respectively.  (*) 

c.l  and  C.2.  As  in  Fig.  2a,  except  that  (m^)  is  changed  from  (0)  to  (0.75). 


Fig.  3.1.  Induced  reactive  factor  -S'(>’)[=  as  a  function  of  (y),  for  a  stiffness 

control  coupling  form  with  =1-^  [a =0.0.],  g  =  =0.  [Sprung-masses.] 

a.  With  6  =  (0.1) 

b.  With  b  =  (2.0) 


Fig.  3.1.  Induced  reactive  factor  5'(;;)[=  as  a  function  of  (y),  for  a  stiffness 

control  coupling  form  with  =1.0  [a  =  0.0.],  g  =  fh^  =^.  [Sprung-masses.] 

c.  With  6  =  (10) 

d.  A  superimposition  of  b  =  (0.01),  (0.1)  and  (0.3),  and  (FOA).  (*) 


Fig.  3.2.  Induced  loss  factor  rj^iy)  [=  ri°{y)],  as  a  function  of  (y),  for  a  stif&iess 
control  coupling  form  with  =1.0  [a=0.0.],  f  =  =  0.  [Sprung-masses.] 

c.  With  6  =  (10) 

d.  A  superimposition  of  6  =  (0.01),  (0.1)  and  (0.3),  and  (FOA).  (*) 


Fig.  3.1.  Induced  reactive  factor  S(y)[=  S^(y)],  as  a  function  of  (y),  for  a  stiffness 
control  coupling  form  with  =1.0  [a =0.0.],  g  =  mc=0.  [Sprung-masses.] 

e.  A  superimposition  of  a,  b,  and  c,  and  (FOA). 

f.  As  in  3  .e  above  except  that  R  =  (7)  and  not  (27). 


Fig.  4.1.  Induced  reactive  factor  5'(>') ,  as  a  function  of  (y),  for  a  stiffness  control 
coupling 

form.  [R  =  27  and  (M,  /M)  =  0. 1 .] 

a.  =  0.75  [a  =  0.25.],  g  -  =0.  [Strong  coupling.] 

b.  =  0.15  [a  =  0.85.],  g  =  fn^  =0.  [Moderate  coupling.] 


Fig.  4.2.  Induced  loss  factor  ri,{y),  as  a  function  of  (y),  for  a  stiffiiess  control 
coupling  form.  [R  =  27  and  (M^  IM)  =  0.1 .] 

a.  =  0.75  [a  =  0.25.],  g  =  =  0.  [Strong  coupling.] 

b.  =  0.15  [a  =  0.85.],  g  =  =  0.  [Moderate  coupling.] 
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Fig.  6.1.  Induced  reactive  factor  S(y),  as  a  function  of  (y),  for  a  mass  control  coupling 

form.  [R  -  27  and  (M,  I M)  =  0.1 .] 

a.  cCc  -S  -  =  1.15.],  =  0.15.  [Moderate  coupling.] 

=  1.03.],  fhc  =  0.03.  [Weak coupling.] 


^c~  S  —  =  1.03.],  =  0.03.  [Weak coupling.] 


Induced  reactive  factor  S(y),  as  a  function  of  (y),  for  mixed  control  coupling 
R  =  27  and  (M^ /M)  =  0.1 .] 
in^  =  0.75  [a  =  1.0.],  g  =  0 .  [Very  strong  coupling.] 

0.53  [a  =  0.47.],  g  =  0.54,  fn^-0.  [Strong  coupling.] 


Fig.  7.2.  Induced  loss  factor  rj^  (y) ,  as  a  function  of  (y),  for  mixed  control 
coupling  forms.  [R  =  27  and  (M^  / M)  =  0. 1 .] 

a.  =  0.75  [a  =  1.0.],  g  =  0.  [Very  strong  coupling.] 

b.  =  0.53  [a  =  0.47.],  g  =  0.54,  =0.  [Strong coupling.] 


Fig.  7.1.  Induced  reactive  factor  *S'(>^) 
forms.  [R  =  27  and  {M^  / M)  =  0.1.] 

c.  =  0.10  [a  =  0.90.],  g  =  0.11, 

d.  =  0.02  [a  =  0.98.],  g  =  0.022, 


,  for  mixed  control  coupling 


[Moderate  coupling.] 
I.  [Weak  coupling.] 


Fig.  7.2.  Induced  loss  factor  rj^^y),  as  a  fimction  of  (y),  for  mixed  control 
coupling  forms.  [R  =  27  and  (M^ ! M)  =  0.1.] 

c.  =  0.10  [a  =  0.90.],  f  =  0.11,  =  0.  [Moderate  coupling.] 

d.  =  0.02  [a  =  0.98.],  g  =  0.022,  =  0.  [Weak  coupling.] 


Fig.  9.  Modification  of  pseudo-statistical  variations: 

c.  In  the  loss  factors  assigned  to  b.  In  the  loss  factors  assigned  to  a.  In  the  resonance  frequency 

individual  satellite  oscillators,  with  individual  satellite  oscillators,  with  distribution  of  satellite  oscillators,  as 

(^>^)  as  depicted  in  Fig.  8c.  [cf  Figs.  {b^)  as  depicted  in  Fig.  8b.  [cf.  Figs.  a  function  A(r),  with  A(r)  as 


Fig.  10.1.  Induced  reactive  factor  S(y)  =  [5*^ (>>)],  as  a  function  of  (>>),  for  a  strong  coupling  defined  by  =  1.0[a  =  0.0.], 
in^  =  g  =  0;  implementing  a  number  of  pseudo-statistical  variations  [R  =  27  and  (M^  /M)  =  0.1 .]  First  order  approximation 
is  superimposed. 
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